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We study a holographic model realizing an "antiferromagnetic" phase in which a global SU{2) 
symmetry representing spin is broken down to a U{1) by the presence of a finite electric charge 
density. This involves the condensation of a neutral scalar field in a charged AdS black hole. We 
observe that the phase transition for both neutral and charged (as in the standard holographic 
superconductor) order parameters can be driven to zero temperature by a tuning of the UV confor- 
mal dimension of the order parameter, resulting in a quantum phase transition of the Berezinskii- 
Kosterlitz-Thouless type. We also characterize the antiferromagnetic phase and an externally forced 
ferromagnetic phase by showing that they contain the expected spin waves with linear and quadratic 
dispersions respectively. 
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I. INTRODUCTION 

Quantum phase transitions naturally occur in strongly 
correlated many-body systems, which often contain com- 
peting interactions and the concomitant competing or- 
ders. When the transition is continuous, or first order 
with weak discontinuities, it gives rise to fluctuations that 
are both quantum and collective [ll-Q- Such quantum 
criticality serves as a mechanism for some of the most 
interesting phenomena in condensed matter physics, es- 
pecially in itinerant electronic systems 0, Among 
these are the breakdown of Fermi liquid theory and the 
emergence of unconventional superconductivity. 

Quantum criticality is traditionally formulated within 
the Landau paradigm of phase transitions. The critical 
theory expresses the fluctuations of the order parameter, 
a coarse-grained classical variable manifesting the break- 
ing of a global symmetry, in d + z dimensions [ij ; here 
d is the spatial dimension, and z the dynamic exponent. 
More recent developments 0,01: however, have pointed 
to new types of quantum critical points. New modes, 
which arc inherently quantum and are beyond order- 
parameter fluctuations, emerge as part of the quantum 
critical excitations. Quantum criticality is hence consid- 
erably richer and more delicate than its thermal classical 
counterpart. In turn, new methods are needed to search 
for, study, and characterize strongly coupled quantum 
critical systems. 
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A. AdS2 and the emergent IR CFT 

A promising new route has come from an unexpected 
source, the AdS/CFT duahty [1], which equates a grav- 
ity theory in a weakly curved {d + l)-dimensional anti- 
de Sitter (AdSd+i) spacetime with a strongly-coupled 
d-dimensional field theory living on its boundary. The 
gravity is classical when the corresponding field the- 
ory takes a large N limit. This maps questions 
about strongly-coupled many-body phenomena to solv- 
able single- or few-body classical problems in a curved 
geometry, often that of a black hole. 

For considering a boundary theory at a finite density, a 
particularly simple gravity setup is a charged black hole 
in AdS which describes a boundary conformal field 
theory at a finite chemical potential for a U{1) charge. A 
rich body of phenomena has been found in this relatively 
simple context including novel transport (see for a re- 
view) , holographic superconductors Il3.[l3l ^ (see |14|-fl6l 
for reviews), and non- Fermi liquids [l7H20j (see also |2l| - 
[29jV In particular, at low energies the system has an 
infrared fixed point described by a (0 -I- l)-dimensional 
CFT [lOl. The CFT has nontrivial scahng behavior only 
in the temporal direction and is represented on the grav- 
ity side by a near horizon geometry AdS 2 x E''^^. This 
realization also offers a unified way, from the gravita- 
tional perspective, to understand both the onset of the 
superconducting instability and the emergence of non- 
Fcrmi liquids. 

One finds that each operator in the UV theory develops 
an anomalous IR scaling dimension S. For example, if a 
scalar operator of charge q and UV dimension A is dual 
to a minimally coupled scalar field in the bulk, its IR 
dimension is given by (see e.g. equation (43) of (20| : in 
this paper wc have set gp to 1) 



1 / A(A-d) 1 
^ 2 ^ V did-1) 2d{d-l) ^ 4 



(1) 



This determines the onset of an instability: if 5 becomes 
complex, the system is unstable^. Due to the term pro- 
portional to —q^ in ([T}, even an operator which is irrel- 
evant in the UV can be unstable in the IR if q is suffi- 
ciently large. This gives rise to interesting examples of 
the phenomenon of dangerously irrelevant operators and 
possibly suggests a novel type of pairing instability driven 
by an IR fixed point. 



^ Technically these should be called charged superfluids, as the 
symmetry that is broken is global and is not gauged within the 
holographic framework. 

^ In the bulk, a complex 5 is a violation of the AdS2 BF bound 
and corresponds to an infinitely oscillating bulk field near the 
horizon and hence an instability, as pointed out in 13(1 . 

^ Note that here a clear understanding of the underlying physical 
mechanism is somewhat hindered by the fact that at this moment 
the gravity model provides only a macroscopic description of the 
condensation. 



A similar formula exists for fermions. Here there are 
no instabilities, but the IR dimension 5 controls the na- 
ture of small excitations around a Fermi surface [lO] and 
in particular whether there exist stable quasi-particles. 
It also naturally yields fermion self-energies that are sin- 
gular only in the temporal direction, properties that are 
characteristic of the electron self-energy in models for 
non-Fermi liquids 32, 33] and the spin self-energy of 
quantum critical heavy fermions @ . 

In this paper we extend these development to model 
itinerant quantum magnetic systems, condensed matter 
systems in which the quantum critical phenomenon is of 
great current interest 
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B. Gravity formulation of a magnetic system 

In the low energy limit of an electronic system, spin- 
orbit couplings become suppressed and spin rotations are 
decoupled from spacetime rotations. Thus spin rotations 
remain a symmetry even though the rotational symme- 
tries may be broken by a lattice or other effects. The low 
energy theory is then characterized by an SU{2) global 
symmetry which describes the spin rotation and a [/(I) 
symmetry describing the charge. This SU{2) symmetry 
is of course only approximate and will be broken at high 
energies. However, if one is only interested in univer- 
sality classes describing only low energy behavior (which 
does not involve spin-orbit couplings), this high energy 
breaking will be irrelevant and one might as well replace 
it by a UV completion in which the dynamics have an 
SU{2) symmetry that is exact at all energies. 

Using the standard AdS/CFT dictionary, the con- 
served currents j^,a = 1,2,3 for SU(2) spin symmetry 
and for U{1) charge should be dual to bulk gauge 
fields making up a S'i7(2)spin x f7(l)charge gauge group in 
the bulk. For simplicity, in our bulk description we will 
consider a theory in which this SU{2) x U (1) symmetry is 
exact to all energies. Modeling spin-orbit couplings and 
understanding how to take into account the relation with 
spacetime symmetries are interesting questions which will 
be left for future study (see also recent discussion in [2^ ) . 

We will be interested in studying the gravity dual of 
an "antiferromagnetic" phase in a continuum limit. In 
such a limit the background value of the spin density is 
zero, but there exists a staggered spin order parameter 
^"^jQ = 1,2,3 which transforms as a triplet under spin 
rotations. Its background value spontaneously breaks 
SU{2) to the [/(I) subgroup corresponding to rotations 
about a single axis. This leads us to introduce a real 
scalar field 0° transforming as a triplet under SU (2)spin 
(and neutral under J7(l)chargc) as the corresponding bulk 
field. A phase with vanishing SU (2) gauge fields but with 
a normalizable 0" ^ in one direction can then be inter- 
preted as an antiferromagnetic (AFM) phase, or a spin 
density wave (SDW) phase in an itinerant-electron con- 
text. A "ferromagnet" would have nonzero SU{2) gauge 
fields, corresponding to a nonzero background spin den- 
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sity in the field theory; we also study this by applying a 
source analogous to an external magnetic field. 

Note that the spirit of this discussion is parallel to 
that of holographic superconductors in that the gravity 
description captures the macroscopic dynamics of the or- 
der parameter and the symmetry breaking pattern, but 
does not explain its microscopic origin. 



C. Quantum phase transitions in holographic 
models of symmetry breaking 

To describe the AFM order and its transition, we pro- 
ceed in a way analogous to the superconducting case. 
The magnetic case considered here involves a neutral or- 
der parameter. Indeed, with the exception of supercon- 
ductivity, most ordering phenomena in condensed mat- 
ter systems involve a neutral order parameter; other ex- 
amples include charge density wave order and Pomer- 
anchuck instability. We are therefore led to consider 
holographic phase transitions involving condensation of 
a neutral scalar field in a finite density system. Here we 
again consider a charged black hole which is dual to a 
boundary conformal theory at a finite chemical potential 
/i for a [/(I) charge. 

Setting q ~ Q in ^ one finds that S becomes complex 

if 



max 



2 



d- Ac] < A < Ac 



Vd 



(2) 



where corresponds to the unitarity bound on a scalar 
operator in d spacetime dimensions. On the gravity 
side this regime is where a scalar field (j) satisfies the 
Breitenlohner-Freedman (BF) bound [s^ of AdS^+i but 
violates the BF bound of the near horizon AdS2 region, 
as was recognized first in [sl] . 

In this paper we construct the phase diagram for the 
condensation of such a neutral scalar field in the T — A 
(equivalently, the T — vn? , see Eq. (18)) plane. For aU 
A in the range ^ there exists a critical temperature Tc 
below which the neutral scalar field condenses. The phase 
transition is second order with mean field exponents. For 
the standard quantization"* in AdS^+i, we find that Tc 
decreases with an increasing A and approaches zero at 

Ae. 

Thus if we allow ourselves to vary A through Ac, we 
find a quantum phase transition, as was previously no- 
ticed in [3^ [3^ in the charged case. Within the context 
of a specific boundary field theory varying the UV confor- 
mal dimension A seems somewhat artificial: however the 
physics here is actually controlled by the effective mass 
of the scalar in the near horizon AdS2 region, and we 



* The story for alternative quantization is more involved, see dis- 
cussion in scc.|Tl]for details. 



later discuss several concrete ways to achieve this. In- 
terestingly, the quantum phase transition at Ac is not 
described by mean field exponents, but is instead of the 
Berezinskii-Kosterlitz-Thouless (BKT) type with an ex- 
ponentially generated scale. More explicitly, for A ~ Ac, 
at T = 0, there exists an IR scale 



Kir ^ /.i exp - 



C 



VAc 



C 



d{d-l) 
"V2A— d 



below which new physics appears (or in other words, the 
condensate becomes significant). From the point of view 
of the near horizon geometry AdS2 , this behavior can be 
understood from the recent discussion in (23l . ISTj which 
argues that for a scalar field in AdS with a mass square 
below the BF bound, such an exponential scale should 
generally be generated.^ An explicit holographic model 
realizing this phenomenon was also recently constructed 
in [s^ and a different quantum phase transition which 
occurs outside the BF boundary (i.e. the BF bound of 
AdS2 is not violated) has been constructed in [39| . 



(0 + l)CFTfl (0+1)CFTa 



AdS: 




UV 



IR 



FIG. 1: A cartoon picture for the flow of the system induced 
by the condensation of a neutral scalar field. The CFTa refers 
to the (O-l-l)-dimensional IR CFT of the uncondensed system, 
described geometrically by an AdS 2 factor with radius i?2. 
When the dimension A of the operator is close to the quantum 
critical value Ac, the system stays near this IR CFT for an 
exponentially long scale, before flowing to the new fixed point, 
(0 -I- 1) CFTb, described by an AdS2 factor with a different 
radius i?2. 



Below the IR scale Ajr we show that the system in 
fact flows to a new IR fixed point which is controlled by 



At the BF bound an IR and UV fixed point, which corresponds 
to standard and alternative quantization respectively, collide and 
move to the complex plane, in analogue with the BKT transition. 
Note that in our case the standard and alternative quantization 
described here refer to those in AdS2. 
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a (0 + l)-dimensional CFT, dual to an AdS2 with a differ- 
ent cosmological constant determined by the condensed 
vacuum of the neutral scalar field. In Fig. [1] we give a 
cartoon picture of this flow. 

The nature of the ordered phase can also be charac- 
terized in terms of its collective modes. The condensate 
breaks the global SU{2) — >■ U{1): we show that the gap- 
less spin waves expected from such a breaking arise nat- 
urally in the gravity description, and that they obey the 
proper dispersion relations. 

The main focus of this current paper is on the con- 
densation of a neutral scalar; however it is clear that one 
can immediately apply the above discussion to conden- 
sation of charged scalar operators which give rise to the 
well-studied holographic superconductors. Equation ([T]) 
implies that dialing A one finds a quantum critical point 
at Ac given by the larger root of (for standard quantiza- 
tion) 



A 



Ac(Ac - d) 

d(rf-l) 2d{d-l) 



1 

4 



(4) 



Again for A close to Ac the system will linger around 
the IR fixed point of the uncondensed phase for an ex- 
ponentially large scale ([3]) before settling into new fixed 
points. In this case the gravity description of the IR fixed 
points have been worked out before for d = 3 in (40l - [43| 
(see also [i^ ) . The stucture of the phase diagram will in 
general depend on the detailed form of the action, and 
non-minimal couplings such as those typically found in 
stringy embeddings (e.g. see [42| - |45j ) will change the re- 
sults. It was found in that for minimally coupled 
scalar actions and potentials similar to ours and for q 
small the system flows to a Lifshitz fixed point with a 
dynamic exponent z ^ -i^, while for larger values of q 
the system flows to another AdS4 with a different cos- 
mological constant. Combined with the discussion above 
for the condensation of a neutral scalar field we thus find 
a unifying picture for the quantum phase transitions for 
both charged and neutral order parameters as presented 
in Fig.m Note that AdS2 can be considered as describing 
a d-dimensional theory with a dynamic exponent z = oo, 
while AdS4 a theory with z = 1. 

Note that in our above discussion we have imagined the 
existence of an "experimental knob" which can be used to 
adjust the UV scaling dimension A of an operator. It is 
important to note that we do this purely for convenience, 
and the nature of our discussion is completely insensitive 
to the precise realization of such a knob. The most useful 
knob will likely depend on the UV geometry into which 
this AdS2 is embedded. For example, if one is studying 
a holographic superconductor an external magnetic field 
will also allow one to tunc the IR scaling dimension . 
A similar UV realization is provided in a D3/D5 brane 
construction in [ssj . where a precisely analogous transi- 
tion is studied. For convenience in the remainder of this 
paper we will simply imagine that we are free to tune 
the bulk mass and thus directly the UV dimension of the 
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FIG. 2: Tuning the UV dimension of the order parameter 
we find quantum phase transitions between a (0-1-1) dimen- 
sional IR CFT corresponding to AdS2 in the unbroken phase 
and various types of symmetry-breaking phases. The type of 
symmetry-broken phase depends on the charge q of the order 
parameter. 



scalar. In sec. IIV Al we give a simple model illustrating 
how this may be achieved. 

A rough plan of the paper is as follows. In the next 
section we describe a holographic phase transition cor- 
responding to the condensation of a neutral scalar order 
parameter in a finite density system in a probe approx- 
imation. In section IIIII we discuss some qualitative fea- 
tures of the effects of backrcaction and the zero temper- 
ature limit of our solution. In section IIVI we discuss the 
quantum phase transition that we find by tuning the UV 
conformal dimension. In section |V] we embed the scalar 
solution discussed in sec.|lI]into an SU{2) system describ- 
ing the spin rotational symmetry. We proceed to study 
the spin waves from spontaneous breaking of the spin 
symmetry in the antiferromagnetic phase. In a probe 
limit we are able to isolate the spin wave excitations di- 
rectly in the bulk and find their dispersion relations. We 
do not construct a spontaneous fcrromagnct; however in 
section IVll we do show that if one considers aligning the 
"spins" with an external magnetic field then techniques 
similar to those in section |V] can be used to find a spin 
wave which has a quadratic dispersion relation, in line 
with field theoretical expectations. We conclude in sec- 
tion EH] with a discussion of further directions. 



II. CONDENSATION OF A NEUTRAL ORDER 
PARAMETER AT A FINITE DENSITY 

In this section we provide a gravity dual description 
of the condensation of a neutral scalar field in a charged 
AdS black hole geometry which describes the onset of a 
real order parameter in the boundary theory at a finite 
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density. We begin our analysis by studying the transition 
at a finite temperature. While our discussion applies to 
any spacetime dimension, for definiteness we will consider 
a (2 + l)-dimensional boundary theory. 



(BF) bound, while still satisfying the AdS4 BF bound, 
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i.e 



- < m'^R'' < - - 



(11) 



A. Setup 

To put the system at a finite density we turn on a 
chemical potential /.j for the C/(l)chargc in the boundary. 
This is described on the gravity side by a charged black 
hole with a nonzero electric field for the corresponding 
U{1) gauge field Bm- The action for B]\i coupled to AdS 
gravity can be written as 



_6_ 



(5) 



with Gmn = QmBm — OnBm and R is the curvature 
radius of AdS. The equations of motion following frorn ([5]) 
are solved by the geometry of a charged black hole [9l.[l0j. 



= gMNUX 



with 



l + 3r; 



Bt 



dx^) + — — 
J 



1 J, 



where we have rescaled the coordinates so that the hori- 
zon is at r = 1 and all coordinates are dimensionless (see 
Appendix lA II for more details). The chemical potential 
and temperature are given by 



^ = ^(1 
47r ^ 



77 is a parameter between and 1, where 77 = 1 cor- 
responds to the extremal black hole with T — Q and 
77 = corresponds to a finite temperature system with 
zero chemical potential (and charge density).^ In the zero 
temperature limit, the near horizon geometry reduces to 
AdS2 X with the curvature radius of the AdS2 region 
related to that of the UV AdSA by 



i?2 = ^ 



R_ 



(9) 



where the lower limit is the BF bound in AdS4, the upper 
limit is the BF bound for the near horizon AdS2 region, 
and we have used © to convert from AdS2 to AdS4 radii. 

Once this condition is met, the scalar will want to con- 
dense near the horizon but will be stable at infinity. The 
condensed solution will involve a nontrivial radial pro- 
file for the scalar; we will see that at low temperatures 
the scalar will probe the extreme values of its potential 
and nonlinearities in the potential will be important. We 
choose to study a nonlinear Mexican hat potential with 
the Lagrangian for x given by 



1 



-\{dx?-V{x) 



with 



m^R^? 



(12) 



(13) 



(6) Here is the effective mass near the point x = 0, and 
should be chosen to satisfy the condition ([TT|) (in partic- 
ular, it is negative). A is a coupling constant and we have 
chosen the constant in (fT2|) so that at x = 0, there is no 
net contribution to the cosmological constant. The pre- 

(7) cise form of the potential in (jlSp is not important for our 
discussion below, provided that it does have a minimum 
and satisfies the condition ((TT|) ^. 

At zero temperature, we expect the scalar to condense 
until the value at the horizon reaches some point near the 
bottom of the Mexican hat, at which point the "effective 
AdS2 mass" will again satisfy the AdS2 BF bound and 
condensation will halt. At finite temperature, we expect 

(8) a phase transition at some temperature Tc below which x 
condenses. Note that at the classical level the coefficient 
of the x^ term is arbitrary, as it can be absorbed into A 
in (fT^ via a rcscaling of x- 



B. Phase diagram 

We now seek the cndpoint of the instability, i.e a non- 
trivial scalar profile for x- We first consider the finite 



As discussed in (20|, |36[, a neutral scalar field x can 
develop an instability if the mass square of the scalar 
violates the near-horizon AdS2 Breitenlohner-Freedman 



Note that since we are considering a conformal theory, only the 
dimensionless ratio ^ is physically relevant, and it is clear that 
as r) is varied from to 1, ^ takes all values from to cxd. 



For general d boundary theory dimensions, equation l llOII be- 
comes 



d2 



< rr?B? < 



d{d- 1) 



(10) 



The condensed phase q = solutions without the stabilizing 
term in the potential were studied in [3^ . |41[ : at nonzero temper- 
ature these are similar to ours, but in the low temperature limit 
the structure described in Section IlIII - which depends critically 
on the existence of a minimum to the potential - is not shared 
by those examples. 
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temperature case and take A to be parametrically large 
so that we can ignore the backreaction of x to the back- 
ground geometry. We wih discuss the backreaction in 
section Hm and there we argue that this approximation is 
good even at zero temperature. 

The equation of motion for x(r) is given by 

^dr (r^fdrx) - X(x' + ^n^R^) = 0. (14) 

We are interested in a solution which is regular at the 
horizon and normalizable at the boundary, which can be 
found numerically. 

We first consider the asymptotic behavior for % near 
the horizon and the boundary. We require the solution 
to be regular at the horizon, i.e. to have the expansion 

x{r) = Xh + x'h{r-l) + --- , r^l (15) 

At finite temperature the factor / in (fT4)) has a first or- 
der zero at the horizon, i.e. /(r) = 47rT(r — 1) + • • • . 
Demanding that be nonsingular then leads to a con- 
dition linking the near- horizon value of x to its derivative: 

x'h = ^Xh{xl + ^^R^) ■ (16) 

A choice of Xh fixes also x'h ^^"^ thus completely specifics 
the solution. 

Near the boundary r oo, the linearized equation 
of (fT4|) gives the standard asymptotic behavior 

X{r) ^ Ar^-^ + Br-^, r ^ oo (17) 

with A given by 



A = ^ + ^m^R^ + ^. (18) 

In the standard quantization A has the interpretation 
of the source^ while B gives the response of the 

order parameter $ dual to x i^i the presence of source 
A. Note that the mass range ([TT|) lies within the range 
— I < m?R^ ^ ~t ^'^^ which an alternative quantization 
exists, in which the roles of A and B are exchanged [48| . 
The conformal dimension of $ (i.e. the "dimension" of 
fluctuations about the point x = 0) is given by A (stan- 
dard quantization) and 3 — A (alternative quantization) 
respectively. 

The condensed phase for standard quantization is char- 
acterized by a normalizable nontrivial solution with A = 
0; and then B gives the expectation value of the order 
parameter <f>. For the alternative quantization we look 
instead for a solution with B = 0, and A then gives the 
expectation value. The task before us now is to pick a 



For example, in section|V]we interpret x ^ corresponding to the 
staggered magnetization, in which case A can be interpreted as 
the staggered magnetic field. 




-0.05 - J 

-oiot ! 

(a) A and _B as a function of the horizon 
value of the scalar field Xh • 

B 

O.lOr 




-O.IOL 

(b) B as a function of A: different points on 
the curve correspond to different values of Xh- 

FIG. 3: The constants A and B determine the behavior of 
the scalar profile asymptotically. This is a representative plot 
where we scan the case m = -2.1 and T = 0.00024 (with 
T/Tc — 0.22) by varying Xh- There is symmetry breaking if 
A = -B 7^ in the normal quantization or if ^ 7^ B = in 
the alternative quantization. 



value for Xh-, and then numerically integrate the radial 
evolution equation to the boundary. As expected for a 
sufficiently low temperature, we find a nontrivial scalar 
hair solution, which means that there will exist a Xh for 
which A or B vanishes, as shown in Figured 

For the standard quantization we find a continuous 
phase transition for all values of falling into the 
range ([TT|) ; there is a nontrivial profile for x for T smaller 
than some temperature Tc and none for T greater. In 
particular, as we increase m^, the critical temperature 
Tc decreases to zero as the upper bound in (fTT|) is ap- 
proached. Precisely at the critical value m^R^ = — |, 
we find a quantum phase transition: the physics in the 
vicinity of this point is discussed in Sec. IIVI The phase 
diagram for standard quantization is plotted in Fig. 2) 

For the alternative quantization, the phase structure 
in the vicinity of the quantum critical point at m? = 
is the same as in the standard quantization. However 
the global structure of the phase diagram is somewhat 
different; in particular the specific value of the mass 
w?R^ — —j^ plays an important role. This value of the 
mass corresponds to the UV scaling dimension for which 
the first nonlinear (i.e. arising from the cj)'^ term in the 
potential) correction to the near-boundary asymptotics 
becomes degenerate with the term proportional to B. For 
inn?R^ < ~fi phase structure is as described above, 
but for rri^R^ ^ ^TE '^^'^ finds a new condensed phase 
in the high temperature regime. In particular, the criti- 
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cal temperature Tc2 appears to increase with m^. These 
solutions appear to be closely related to the thermody- 
namically unstable scalar hair solutions constructed in 
[49I . [50| , which studied uncharged black holes and thus 
correspond to the high-temperature limit of our construc- 
tion. As these new phases appear to involve UV physics 
and are not related to the low temperature quantum- 
critical behavior that is the focus of this work, we defer 
an in-depth study of these phases and the critical value 
w?B? = — yI to later work. 



C. Critical exponents 

A continuous phase transition can be characterized by 
various critical exponents: in the normal quantization 
with A the source and B ^ ((f>) the expectation value we 
have the following behavior close to Tci 

1. The expectation value ($) cx (Tc — T)*^ with mean 
field value P = \- 

2. The specific heat: C cx [T^ - T|"" with mean field 
value a = 0. 

3. Zero field susceptibility: ^^-\a=o oc jT — Tc\~^ 
with mean field value 7 = 1. 

4. Precisely at T = Tc, (<i>) ~ with mean field 
value (5 = 3 



FIG. 5: Top: Plot for exponent /3, defined as B ~ {Tc - Tf. 
P = 0.49±0.03 from numerical fit, compared with /3mean fleW = 
i. Bottom: Plot for exponent 5, defined as B ~ at Tc. 
5 = 3.03±0.05 from numerical fit, compared with i5nioan field = 
3. 



The existence of mean field exponents in this classical 
gravity analysis is not surprising^'^ and can be translated 
into the statement that A and B arc analytic functions 
of the horizon field Xh- This follows from the fact that x 
is regular at the horizon and the black hole geometry is 
smooth; we will find a rather different situation at T = 
in the next section. 

We now sketch the relevant arguments for finite-T. As 
argued above, for small Xh we can expand A as 



A^ai{T)xh + a2{T)xl + ' 



(19) 



and similarly for B. Note that both A and B must both 
be odd functions of Xh, due to the x ^ symmetry. At 
Tc, a new zero of A should be generated at = and 
move to a finite value of x/i as T is further lowered. ^"'^ 
Thus near Tc we should have ai(T) = a(T — Tc) and 
a2 (T) = b + ■ ■ ■ with a and h having the same sign for 
T > Tc- For T < Tc, then A has a zero for Xh taking a 
value XA given by 



XA 



a2(T) 



(X (Tc-T)^ 



(20) 



Using xa a-s the initial value at the horizon gives the 
solution for the condensed phase. Assuming that B is 



There are also other exponents associated with correla- 
tion functions at finite spatial or time separation which 
we will leave for future study. For our phase transition we 
find that all of the above exponents are precisely those of 
mean field theory. Sec fig.[5]for numerical fits of exponent 
/3 and S. 



Mean field exponents also appear in the phase transition associ- 
ated with a charged scalar field in the holographic superconduc- 
tor story [Ulil. 

In fact given x ~^ ~X symmetry, a pair of zeros are generated 
and move to opposite directions. 
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still linear in Xh near Xh = 0, we thus find that 

B(xA)(xXA(x(T,-r)^ (21) 

which gives us the critical exponent /3 — ^. Turning now 
to the susceptibility, we find 



dB 
dA 



dB dxh 



A=0 



dxh dA 



(22) 



leading to 7 = 1. Furthermore, precisely at T = Tc, 
A ^ ax\ and B stays linear near Xh = 0. Thus we find 
that at T = T^, 



B{xh) cx x/i oc ^3 



(23) 



which gives the critical exponent 5 = 3. One can also 
compute the free energy using holographic renormaliza- 
tion and indeed find mean field exponent F oc —{Tc~T)^ 
for T < Tc- In particular, the analytic expansion (jl9p 
for A and B in terms of Xh guarantees that F has the 
Landau- Ginsburg form. 

Note that the gravity analysis can be extended to all 
spacetime dimensions with the same mean field scalings. 
Presumably this has to do with the fact we are working 
in the large N limit, which suppresses fluctuations. This 
is also consistent with picture obtained in (20j in which 
the scalar instability essentially follows from a RPA type 
analysis. We expect that a 1/N computation involving 
quantum corrections in the bulk will reveal corrections 
to these mean-field exponents. 



III. BACKREACTION AND THE ZERO 
TEMPERATURE LIMIT 



hi the usual studies of holographic superconductors, 
the scalar is charged under a U{1) and so the usual 
probe approximation involves treating the scalar and 
U{1) gauge field as negligible perturbations to the back- 
ground metric. In these models lowering the temperature 
essentially means increasing the ratio /i/T; thus at suf- 
ficiently low temperatures the gauge field (and thus also 
the condensed scalar) will necessarily backreact strongly 
on the geometry, causing a breakdown of the probe ap- 
proximation. The zero temperature limit of the backre- 
acted geometries have been constructed jiol - lisj and typi- 
cally depend on the details of the couplings and charge of 
the scalar. One recurring theme is that at zero tempera- 
ture the charges that was previously carried by the black 
hole is now completely sucked out of the hole and into its 
scalar hair. Once the black hole horizon is relieved of the 
burden of carrying a large charge, it is usually replaced 
by a degenerate horizon with vanishing entropy. 

In our model the situation is different. We always in- 
clude the backreaction of the gauge field on the metric, 
as we start from the beginning with the charged black 
hole solution; it is consistent to solve for a scalar profile 



on this fixed background only because the scalar is un- 
charged and so its contribution to the backreaction can 
be cleanly suppressed by taking A large. 

In this section we will discuss the backreacted solu- 
tion in the IR at zero temperature and argue that even if 
backreaction is included its effects are rather benign and 
do not change any qualitative conclusions. This is essen- 
tially because all of the charge must stay in the black 
hole itself, greatly constraining its near-horizon form. 

First, we note that Gauss's law states that 



^gg"g''drAt 







;V5 



"g^'drAt 



±t = const . 

(24) 

Thus if the electric field is nonsingular at the horizon, 
gxx must be finite there: this is simply saying that the 
at the horizon cannot degenerate as it has a nonzero 
electric field flux through it. We then expect that the 
near- horizon geometry factorizes into the form A^2 x R^, 
where M2 is some 2d manifold involving {t, r). Now con- 
sider the trace of the Einstein equation arising from the 
variation of ([5|) plus p2|) : the C/(l) fleld strength does 
not contribute to this equation as it is classically scale 
invariant, and we flnd: 

7^ + ^ = |^[(Vx)2+4y(x)], (25) 

where TZ is the Ricci scalar of the geometry and comes 
purely from the M2 factor. Now let us further assume 
that in the near-horizon region x asymptotes to some 
constant value Xh (we will show this to be consistent 
shortly). We then flnd that A^2 has constant negative 
curvature, and thus must be AdS2^^ with radius R2 sat- 
isfying 



1 
If 



1 



v{x) 



(26) 



where R2 is the curvature radius of the original AdS2 near 
horizon geometry of the extremal charged black hole. 
Note that since V{xh) < 0, i?2 < R2- Thus we see 
that the backreacted IR geometry is very similar to the 
unperturbed geometry, except that its AdS2 factor has 
a radius that is corrected by the presence of the near- 
horizon scalar potential. 

Let us now study the scalar equation of motion 
near the backreacted AdS2 horizon: 



\dAr'fdrX) = R'^ 
dx 



(27) 



where / is the warp factor for the backreacted geometry 
and at zero temperature has a double zero at the horizon 



More precisely, it could also be an AdS2 black hole; the true zero 
temperature solution corresponds however to pure AdS2. 
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rp; expanding near the horizon we see that regularity at 
the horizon requires 

^(X(r = ro)) = 0. (28) 

Thus we see that at the horizon x will sit at the bot- 
tom of its potential. To understand how this AdS2 re- 
gion matches onto the asymptotic geometry, we expand 
X — Xh + 5{r) where Xh is the bottom of the potential 
V{x)- Now working in the AdS2 region and linearizing 
([27| near Xh we find that 5 obeys the standard ^^52 
wave equation: 

dr{{r - r^fdrS) - RlV"{xh)S = 0, (29) 

whose solutions are 

S = a{r^ro)-^+''+P{r-ro)-^-' (30) 

with 

^ = \ll+RlV"{xh), (31) 

Now we would like x to approach Xh as we approach 
the horizon; this means that we should not take the j3 
solution above, as it invariably blows up as r ^ ro. Note 
that for the a solution to also not blow up, we need 
+ > and thus V"{xh) > 0. To stabilize the 
AdS2 region we must truly be sitting at a minimum of 
the potential at the horizon. In this case the solution 

x{r)=Xh+a{r-ro)-^+-' (32) 

can be interpreted as an irrelevant deformation of the 
new AdS2 IR CFT that we are flowing to"^^. The value 
of the coefficient a is not fixed at this linearized level and 
must be determined by matching to the UV solution. 

We have not solved for the fully backreacted geometry 
but have used this method to find a normalizable scalar 
solution on the T = charged black hole geometry using 
the above exponents. Our results match smoothly onto 
the T — >■ limit of the profiles calculated using the finite 
temperature matching procedure discussed in Section HIl 
When A is large we do not expect the inclusion of the 
backreaction to qualitatively change any of these results 
given the change of the cosmological constant is small. 

It is instructive to compare this to the situation for 
charged holographic superconductors [1^. In this case 
for large charge it is found that the IR geometry flows to 
an AdSi, while for sufficiently small charge it flows to a 
Lifshitz geometry, with exponent z satisfying 




Similar considerations are used in [40l | to determine when an 
emergent AdS4 can exist. 



at very small q (see Eq.(81) in [131). We see that our 
AdS2 solution corresponds to z = oo at 9 = 0; increasing 
the charge we have a Lifshitz solution with finite z, and 
finally increasing further we find z = I for AdS^. 

IV. A QUANTUM PHASE TRANSITION FROM 
CLASSICAL GRAVITY 

We recall from the discussion in sec. Ill Bl that when 
m?^B? = — |, the critical temperature approaches zero: 
thus we should obtain a "quantum phase transition" at 
zero temperature as is varied from above to below. 
In this section we consider the behavior near the quantum 
critical point from the condensed side. 

The accessibility of this quantum critical point depends 
on having the ability to tune the IR dimension of the 
field X in some way. The simplest possible way is by 
directly tuning the mass, and before proceeding we pause 
briefly to explain how this may be possible in a simple 
bulk realization. A similar mechanism is discussed for 
uncharged black holes in AdS in [s^ . 

A. Tuning across the quantum critical point 

Imagine that the scalar x contains a coupling to an- 
other scalar field V', 

S^^ j d''+'xV^F{^)x^ (34) 

where Fl-tp) is some function whose detailed form will 
not be important for us. Now consider turning on a con- 
stant source h for ijj at the boundary, i.e. we require the 
boundary value of ip to approach h. In the boundary 
theory language this corresponds to adding a term 

J d'^xhO{x) (35) 

to the action with O the boundary operator dual to -0. 
The source induces a nontrivial bulk solution for which 
in turn contributes to the action of x as a mass term, 
shifting the effective m^. However the new effective mass 
will typically be radially dependent; we require it to be 
nonzero at the AdS2 horizon. If we write the IR dimen- 
sion ll]) of ?/) as 

5^ = i + (36) 

then the asymptotics of ^ in the near horizon AdS2 region 
are 

0^ (r-ro)-5±'^^ . (37) 

For -0 to be regular at the horizon we need to choose 
the -|- sign in the exponent. So if z^^ = ^ then ij) ap- 
proaches a constant at the horizon, and turning on a 
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source h for tp will change the effective IR mass of the 
field X- From ([1]), i/^ is i if V' is dual to an operator 
with UV dimension A = d, i.e. it is masslcss in the bulk. 
Thus whenever the UV CFT has such an operator with a 
nontrivial OPE with the order parameter operator $, we 
expect to be able to tune the system through the quan- 
tum critical point by varying h in ()35p . i.e. there exists 
a critical value he of h at which we expect a quantum 
phase transition. It would be interesting to understand 
this mechanism further. Also note that for type II theo- 
ries in an asymptotic AdS geometry, a natural candidate 
for ip is the dilaton. 



B. BKT scaling behavior and Efimov states 



For 771^ < m^, the BF bound of the near horizon region 
AdS2 region is violated. It has been argued in [s^l that 
for a general AdS gravity dual with mass slightly below 
the BF bound conformality is lost and an IR scale A/j^ 
is generated exponentially as 



A 



/fl ~ /i exp 



(38) 



where represents some UV scale that in our case is the 
chemical potential. This exponential behavior is char- 
acteristic of the Berezinskii-Kosterlitz-Thouless (BKT) 
phase transition. One expects that this IR scale controls 
the physics near the quantum phase transition. In par- 
ticular, the critical temperature Tc and the value of the 
condensate ($) should be related to this scale. Here we 
present a reformulation of the arguments of (STj that 
demonstrates the behavior of Tc and ($) explicitly. 

We parameterize the AdS 2 region of the uncondensed 
geometry as 



(39) 



The linearized equation for the scalar about the point 
X = is then 



dz 



Rl{m^ - ml) - 1/4 



(40) 



where we have assumed time dependence e~*"*. This is 
essentially a Schrodinger equation with "energy" oj^: the 
existence of a negative energy bound state indicates a 
mode growing exponentially in time and hence an insta- 
bility, and it is a well-known fact that if z can take values 
from to 00 then such negative-energy bound states exist 
when < m^, leading to the BF bound. 

In our problem, however, z docs not take values on 
the whole half-line. At some small UV value zjjv we 
should match onto the asymptotic UV geometry. We will 
now show that if we also choose an appropriate infrared 
cutoff zjji then even when is below the BF bound 
(pni) will have no bound states. To understand this, let us 



assume for simplicity that x satisfies Dirichlet boundary 
conditions xi^iR.uv) = and consider the zero energy 
a; = solutions to (j40|) . We find that the solution is 
oscillatory (see also discussion in [23j ) 



Xiz) 



1 z sm 



loe 



(41) 



and importantly, the zero energy solution satisfies the 
boundary condition only when 



lOE 



zuv 



(42) 



Essentially we must fit a single half-period of the oscil- 
latory wave function inside. This wave function has no 
nodes; thus if zib. satisfies this condition, then the lowest 
energy state has zero energy and there is no instability. 
Decreasing the distance between zi^ and zjjv will only 
increase the energy of the ground state. On the other 
hand, if we increase this distance-if z/_r is too high-then 
the ground state will have negative energy, indicating an 
instability. 

Now in our problem zm can be provided by a small 
finite temperature, which ends the geometry with a hori- 
zon at some large value Zh- We thus find that if Zh > zm 
as defined above, there will be an instability that can be 
resolved only by condensation of the scalar. The critical 
temperature Tc is thus given by 



Zh 



^ exp 



(43) 



where ^ ^ exp[z[/y] is some UV scale. We have been 
able to confirm this numerically, including the prefactor 
in the exponent (see Fig. [6]). 



Tc exp ^— V6 TT 



40- 



30- 



20- 



10- 



-2.2 -2.1 
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-1.9 
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-1.6 



-1.5 



m 



FIG. 6: The exponential dependence of Tc as a function of 
\/m'i — m-^i?2. Here we liave compensated for the expected 
theoretical behavior of Tc and plotted the resulting function 



versus m ; we see that as m — >■ 
constant, verifying our prediction. 



the values approach a 



Note that more general boundary conditions simply result in 
slight changes in the value of zjn and zi/y 
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What if T < Tc? Now the IR cutoff cannot be pro- 
vided by the horizon, and must instead be provided by 
nonUncaritics associated with the scalar potential, which 
will become important near z/^. We then expect that 
the condensate ($) should itself exhibit similar exponen- 
tial scaling as (|43|). as was found explicitly in a similar 
context in (38| . We present here a simple way to under- 
stand this result. Note that for rri^ only slightly below 
TO^ the IR theory is basically the original conformal IR 
CFT for an exponential hierarchy of scales. In this IR 
CFT, X is unstable and has dimension given by 



1 



+ i 



1 



(44) 



Since x has dimension i, 
scaling transformation as 



it should transform under a 



x(Az) = \^x{z) 



(45) 



We expect that ($) should be simply related to the value 
of x{z) the UV matching value zjjy, as essentially only 
linear UV evolution will be required to relate them, and 
thus 



X{zuv) 



(46) 



Similarly, the value of </> in the deep IR will be determined 
by nonlinearities, which we expect to result in 



Xizin) ^ 0(1) 
We thus find from (HSI) that 



x{zir) 



ZIR 

zuv' 



x{zuv) 



(47) 



(48) 



which results upon evolution through the UV region to 
the AdS4 boundary in 



Zuv 
ziit 



Am cxp 



(49) 



Note that the exponent for ($) is half that obtained for 
Tc- In the intermediate conformal regime x scales as a 
dimension 1/2 operator and the temperature scales as di- 
mension 1 : this is the origin of the factor of two difference 
in the exponents above. Eventually in the far IR nonlin- 
earities become important and a scale is generated, but 
the theory is not gapped; provided that the scalar poten- 
tial has a minimum, as described in Section [IIII it instead 
flows to a different AdS2 fixed point, and the IR scale 
manifests itself only as an irrelevant perturbation ([32]) 
along which we flow to this new fixed point. 

It is possible to perform a more explicit calculation 
of (Unili following the techniques used in [s^ where sim- 
ilar scaling is obtained in an AdS2 that arises from a 
D3/D5 brane construction. From such a treatment it is 
clear that there should exist an infinite tower of "Efi- 
mov states" that correspond to allowing oscillatory so- 
lutions such as (|4ip to move through more periods be- 
fore matching to the IR solution. These have ($) ~ 



. with n a positive integer . We have been 

able to find the first two of these states numerically (see 
fig. [7]); as explained in |38j the relevant wavefunctions 
have more nodes with increasing n and so the n = 1 
state is energetically favored. 




FIG. 7; Efimov states: for comparison we plot expected 
results from the analytic arguments described in the text. 
Agreement can be seen at this level but numerical difficulties 
prevent us from probing this region more carefully. 



To conclude this subsection we note that turning on 
any finite temperature the phase transition becomes that 
of the mean field, since the physics depended in a smooth 
way on the horizon value of the scalar Xh- At zero tem- 
perature since the horizon is degenerate we expect the 
boundary values of x to depend on the initial value a at 
the horizon in equation p2p . It would be interesting to 
understand this better. 



Quantum critical points for holographic 
superconductors 



Finally, note that the scaling behavior (|43)) and (j49|) 
also applies to the condensate of a charged scalar, which 
results in the well-studied holographic superconductor.^^ 
In such a case the coupling to the background electric 
field is important; also one now has the possibility of sup- 
plying an external magnetic field H for C/(l)charge^^ j47l |. 
The effect of the electric field on the conformal dimension 
is shown in ([1]) , and the generalization of this formula to 
include the effect of the magnetic field is (see Appendix lAl 



Similar exponential dependence of the critical temperature in 
such a ease has been observed by Faulkner and Roberts. 
Note that this magnetic field is an external field strength for 
(7(l)charge and has absolutely nothing to do with the "magnetic" 
field associated with the antiferromagnetic ordering discussed 
later in this paper, which would correspond to a chemical po- 
tential for the SU{2) gauge field Af. 
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FIG. 8: Contour plot of the critical magnetic field in the m^, 
q plane. The color indicates the value of He, with more purple 
indicating a higher field. In region A the IR dimension even 
with H = Q 18 real and there is no condensate. In regions B 
and C there exists a finite magnetic field He above which the 
condensate is destroyed. However H^ diverges as we move 
to the left and is infinite on the boundary between regions 
C and D. In region D even an infinitely strong magnetic 
field will not stop the condensation. This is understandable: 
everywhere to the the left of the line = the scalar 
has a sufficiently negative mass that it would have condensed 
even if it was uncharged and so it is perhaps not unexpected 
that a magnetic field cannot halt this condensation. In region 
C again the scalar would have condensed even if it was un- 
charged, but here the charge is high enough that a sufficiently 
strong magnetic field can stop the condensation. 



for a derivation) 




i6\bq\ - g2 



Vi + 1252 



- + i (50) 



where b = —5- is dimcnsionless ratio between the bound- 
Ms 

ary magnetic field H and chemical potential^'^. Again the 
quantum phase transition will happen when S becomes 
complex. It is clear from here that there is a different 
critical mass than in the neutral case; alternatively, 
one can now imagine tuning the magnetic field through 
a critical value He, which is again found by setting the 
expression inside the square root of ()50|) to zero. For 
H > He there will be no condensate, and for H sUghtly 



^"^ fiB is the dimensional version of the dimcnsionless chemical po- 
tential /i we were using earlier, see Appendix lA II for further 
explanation. 



less than He one will find similar non-analytic behavior 
as above as a function of the deviation of H from He- 
The explicit expression for He is rather complicated and 
is given in (|A28|) . 



We plot the critical values He in the q 



plane in 



Fig. [H Note in particular that the expression in ([50 
containing b saturates at a finite value 
Thus if 



2^3 



^/3|g| < -- 



(51) 



no matter how large the magnetic field is, a condensate 
cannot be prevented. This is region D in Fig. [51 This 
is understandable: everywhere in this region the scalar 
mass is below the neutral AdS2 BF bound m?B? = — | 
and so it would have condensed even if it was uncharged. 
As the electric field is in some sense not responsible for 
the instability, it may make sense that a magnetic field 
cannot halt it. It will be interesting to find field theoret- 
ical models with this feature. 



ANTIFERROMAGNETISM AND SPIN 
WAVES 



We have understood in detail the mechanism by which 
a neutral scalar field can condense in a finite-density ge- 
ometry. We would now like to employ this new under- 
standing to holographically model symmetry breaking 
with such a neutral order parameter. One immediate 
example is antiferromagnetic order: as explained in the 
introduction, in the low-energy limit of interest to us spin 
rotations may be thought of as a global SU{2) symmetry. 
Antiferromagnetic ordering then corresponds to a spon- 
taneous breaking SU{2) — >■ U{1), where the unbroken 
direction corresponds to rotations about a single axis. In 
such a system the background value of the spin density 
remains zero, but there is a neutral order parameter 
that corresponds to the staggered magnetism. 

In this section we construct the bulk dual of the op- 
erator $°: we embed the neutral scalar field x discussed 
in sec. [IimVl into part of an SU{2) triplet scalar field 
charged under the SU (2) bulk gauge field corresponding 
to the global spin SU{2) symmetry of a boundary the- 
ory. The phase transition discussed earlier then becomes 
a transition to an antiferromagnetic (AFM) phase. 

Note that while we use the word "antiferromagnetic" 
there is no sense in which the microscopic degrees of free- 
dom of our system consist of spins that arc anti-aligned 
on a bipartite lattice. We use the term to describe only 
the symmetry breaking pattern described above, realized 
in a manner that manifestly requires a finite density for 
symmetry breaking. We do feel that if indeed gravity du- 
als could be constructed top-down for such systems they 
would likely contain ingredients similar to those in our 
description. 
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A. Embedding of x 

More explicitly, we can consider the following action 

S = Sgiav + >5'matter (^2) 

with 



1 6 



2^r^^r+^) (53) 



where R is the curvature radius of AdS4. The relevant 
part of the bulk Lagrangian for matter fields is then given 
by 



C 



matter 



T^a ipMNa ri r<MN 

1 /I 



A V 2 



with 



and 



^NiN — dAiA^ — On All + e°^^^A/^5^, 
Gmn — QmBn — OnBm, 



(54) 



(55) 



Dm^ = + . (56) 

We will take the potential V to have the double well form 
of (HI 



F(0") = JmV-<? + 



(57) 



We again put the system at a finite density by turning 
on a chemical potential for the C/(l)chargc, with a back- 
ground metric given by ©. We note that in this back- 
ground all SU{2) gauge fields are inactive. This high- 
lights an important difference between this setup and the 
usual U{\) holographic superconductor; in the Abelian 
case there is a background chemical potential /i for the 
[/(I) charge; this does not break the U{\) but does in- 
teract via the bulk equations of motion with the charged 
scalar order parameter, causing it to condense for a suit- 
able choice of mass and charge. This is distinctly different 
from the SU (2) case studied here; specifying a chemical 
potential for the SU(2) would involve picking a direction 
/i" in the SU (2) space, corresponding to an explicit break- 
ing of the symmetry. This is analogous to applying an 
external magnetic field. We study this explicit breaking 
in sec. lVII and for now focus on the spontaneous breaking 
of SU{2). 

When falls in the range ([TT|). the background © 
becomes unstable towards condensation of the scalar 0", 
which will then spontaneously break the SU (2) to a U{1) 
subgroup. More explicitly, consider the following ansatz 



= 0,0, 



R 







(58) 



which can be readily checked by examining the equations 
of motion as being self-consistent, if we ignore the back- 
reaction to the background geometry. We note in partic- 
ular that one can consistently set the SU (2) gauge field 
to 0. This is because of the non-Abelian nature of the 
interactions of the gauge field and scalar, e.g. f°'^'^4)bAc 
etc. which clearly vanish if all objects point only in one 
direction in the SU{2) space. 

Plugging ((58)) into the equations of motion following 
from ()54p we precisely find ()14p and our previous discus- 
sion of the phase transition in sec. ITll4IVI can be taken 
over completely. Note that if we consider a finite A and 
turn on backreaction, the profile for the U{1) gauge field 
in ([5]) and the background metric will be modified, but 
one can still set the SU (2) gauge field to zero. 



B. Spin waves 

To characterize the ordered phase, we study in this 
section perturbations around the symmetry breaking so- 
lution ()58|) with a nonzero order parameter 0° but zero 
background gauge field Afj. We show that the system 
has two linearly dispersing gaplcss modes, correspond- 
ing in the field theory to the two Goldstonc modes aris- 
ing from the spontaneous breaking of the global SU{2) 
symmetry, and that their velocity is given by the expres- 
sion expected from the standard field theory for a quan- 
tum antiferromagnet. From the bulk point of view, low- 
frequency rotations of the order parameter in the broken 
symmetry directions will source bulk gauge field fluctua- 
tions, which we will find to be normalizable at the AdS 
boundary if they obey a specific dispersion^^. 

Our analysis for the remainder of this section will not 
depend at all on the details of the metric or scalar profile 
(except that it is normalizable) . We will work in a general 
boundary spacetime with a bulk black brane metric given 
by 

ds^ = gMNdx'^^ dx^ ~ gttdt^ + Qrrdr^ + QxxdxP' (59) 

We work at finite temperature with gu (grr) having a 
first order zero (pole) at the horizon r = tq. 

Let us begin by understanding why from the gravity 
point of view there should exist a gapless mode. Consider 
a global SU{2) gauge rotation of the background order 
parameter 4>Q{r), at lu = and fc = 0. 



<5C(0 = era,4{r) 



(60) 



where r' , i = 1 , 2 is a generator along one of the broken 
directions so that S(f>'^ is nonzero and e is a constant. This 
small perturbation will obviously be a solution to the 



Methods similar to those used here may bo used to obtain an- 
alytic control over the hydrodynamics of the U(l) holographic 
superconductor |55|1 . 
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bulk gravity equations of motion; it is in fact a normal- 
izable solution at the AdS boundary, precisely because 
(pQ is normalizablc. This appears somewhat trivial — but 
note that this is exactly a gapless mode, as it is a nor- 
malizable solution that exists in the limit w — >■ 0. Now 
consider a local SU (2) gauge rotation with the gauge pa- 
rameter e having a small frequency and momentum in 
the x-direction, i.e. e{t,x) oc e~"^*+''^^. The perturba- 



tion to the scalar takes the same form (|60|) . but now we 
have extra perturbations to the gauge fields 



5AI 



ike 



(61) 



These arc not normalizablc at infinity, and thus this pure 
gauge transform is no longer a normalizablc solution. 
Thus at any nonzero u and k to find a normalizablc mode 
we must move off of the pure gauge solution and actually 
solve the dynamical equations of motion. The existence 
of the the global solution ((60|) guarantees that we as we 
take w, fc — > we will find a gapless solution, but the dy- 
namical bulk equations of motion will show that it will 
be normalizablc only when a certain dispersion relation 
ijj{k) is satisfied. 



The equation of motion for tt is then 

dM {x^V^g^'^idNTT - An)) = (65) 

It is clear that in the limit w. A: — > 0, a constant pion 
profile 7r(r) = ttq and vanishing gauge field A ~ provide 
a solution to the equations of motion. This is the gauge 
transform discussed earlier. 

Let us now briefly detour to discuss the asymptotic 
behavior of the pion field. Carrying out the standard 
analysis and using the fact that the background solution 
is normalizablc^'': xi^") ^ '""^ for large r, we find that 
as r — > oo. 



7r(r) K B + A 



,-d+2A 



r — >■ oo 



(66) 



To understand this it is useful to remember that fiuc- 
tuations in the original field S(f)(r) are related to tt as 
S(j){r) = 7r(r)V*0o(?') ~ 7r(r)r~^. Comparing this to the 
asymptotic behavior of the scalar ([TT]) we see that the 
constant piece B of tt at infinity is actually normaliz- 
able, while the piece A is not normalizablc and should be 
viewed as the source for the pion field. 



1. Pions in the bulk 

The previous discussion has convinced us that the exis- 
tence of the gapless mode is intimately related to our abil- 
ity to perform global rotations on the order parameter, 
i.e. to bulk Goldstone modes. Let us thus parametrize 
fluctuations of 4> in terms of bulk Goldstone fields (or 
"pions") 7r'(r, t,a:) as follows: 



4>{r,x) = exp{iTr^(r,t,x)T^)(j)o(r) 



(62) 



2. Yang-Mills equations 

Fluctuations of tt will excite the non-Abclian gauge 
field A'^j(r,t, x), whose linearized equations are 
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V M-t 



xV^(ap7r-^p) 







(67) 



Again we have suppressed the index i which labels the 
broken generators, since at the linearized level different 
directions decouple. We have absorbed various factors 
into a rescaled gauge coupling 



where (f>o{r) is the background solution and as before r* 
are the broken symmetry generators. We now work out 
the bulk quadratic action of the tt* to be^^ 



1 



9 a 



(68) 



2k2 / 2A 



Sin) 

ZK- I ZA 

(63) 

Here hij can be viewed as a metric on the Goldstone 
boson space, given by 



%M = l4'Ur){n,T,}Mr) - j^S^j (64) 

Thus the two Goldstone modes decouple. Below we will 
focus on one of them and drop the index i for notational 
simplicity. 



which controls the strength of the effects of the scalar sec- 
tor on the gauge field. We could of course always choose a 
(Sa/tt — ^A/)(97v7r^— ^i/)/iixinitary gauge to get rid of the Goldstone boson tt — but 

to make the connection to boundary Goldstone modes 
more obvious wc will not do this and instead choose the 
gauge Ar = 0. It is convenient as in [s^ to work with 
the momenta conjugate to the bulk gauge field A^, 
defined as^^ 



1 



~9\ 



-gg^g^^'^drA, 



(69) 



In writing this expression we have assumed that the background 
gauge field is zero; thus we have neglected intrinsically non- 
Abelian terms of the form /°'"^7r(,j4c which all arise at higher 
order. 



We restrict ourselves to the ordinary quantization for this section. 
■^'^ It is amusing to note that computing the momentum conjugate 
to TV precisely extracts out the source in this case, in exactly the 
same way that it extracts out the vev in the more familiar case 
of a standard massless scalar. 

For convenience of discussion here we have chosen a diff'eront 
normalization for the currents below. 
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The value of the bulk field J'' at the AdS boundary r — )• 
oo is equal to the expectation value of the field theory 
current (J''^)qft- 

The N = r component of the Maxwell equations (p7)) 
can now be written 

d^J^ = -V^.9'-'^x'(r)9.7r . (70) 

If the symmetry is unbroken then x = and this is noth- 
ing but the conservation of current. Let us now consider 
evaluating this expression at the AdS boundary. Com- 
paring ([70]) to the asymptotic behavior of tt in ([66|) we 
see that the entire right-hand side of this expression is 
proportional (with no extra factors of r) to the coeffi- 
cient A in the near-boundary expansion of tt, i.e. to the 
source for the Goldstone boson field. This is simply the 
usual field theory Ward identity, which says that in the 
absence of a Goldstone source the field theory current is 
conserved. 

The other nontrivial equations from (|67|) are those in 
t and X directions with nonzero A^ir, t, x) and At(r, t, x) 
only, given by 

-dr,P + ^V^F-' + x/^g"x'(9t7r - At) = (71) 
9 a 

-drr + ^\/tF*^ + V^g'^'^x^d^n -A,) = 0. (72) 
9a 

Equations (|55|) and (|7n)) - ([7^ are the full set of equations 
for this system. 

3. Boundary Goldstone modes and their spin wave velocity 

We will now solve the above equations in the hy- 
drodynamic limit to show explicitly the existence of a 
Goldstone mode, i.e. a normalizable solution to ((65)) 
and ([7n|) - ([7^ that is infalling (or regular) at the hori- 
zon and exists in the limit of small a; or fc. We work in 
Fourier space 

TT 7r(r)e-''^*+*'=^ ^ = At,^(r)e-*"*+''^^ (73) 

Recall that at precisely frequency and momentum we 
already know a solution to these equations: a constant 
pion profile 7r(r) = ttq and A = 0. We now expand around 
this solution in powers of oj and k. To lowest order we 
simply solve the forced equations for A given by (|7T|) and 
(j72p with the forcing term provided by the zeroth-order 
solution for tt, meaning that 

At{r) ^ OiuTTo), A.,{r) ^ 0{^^k). (74) 

The first correction to 7r(r) enters through ([55]) and in- 
volves one extra field theory derivative, and thus we see 
that 

7r(r) = TTo -I- 7ri(r) 7ri(r) ~ ©(ttqw^, ttqA;^) (75) 

The solution we want for A is both infalling at the black 
hole horizon and normalizable at the boundary; such a 



solution is nontrivial because of the forcing term^'^ and 
is given by: 

At = -iujTTo{l - at{r)) , A^ = ikTTii{l ~ a^ir)) (76) 

where at and Gx are defined to be infalling at the horizon 
and satisfy the homogenous part (and zero frequency) 
part of ((7T|) and ([7^ : 

^dr{V^9^'g"drat) - V~99'\^at = (77) 
9 a 

^driV^gg^g^^drax) - V^gg^^x^a, = (78) 
9 a 

To ensure that At and A,j. in (|76p are normalizable we 
also require that 

at(oo) = 1, ax{(yo) = 1 . (79) 

The corresponding equations for the pion profile are 

TTi = w27roC*(r) + k^iroC-'ir), (80) 

where the dynamical equations for C* and can be 
written from (|65p as expressions for the radial indepen- 
dence of two objects a* and a^: 

dra'= = 0, = -yf^x^g^'^'drC'' + ■^yf^g''''g''''drax 

9a 

(81) 

dra' = 0, a* EE sf^x^g^drC - ^^g"g''drat 

9 a 

(82) 

The precise form of these equations turns out to not be 
important: the crucial fact is that C*'^ are essentially 
driven by a^^^, and so we can find solutions to them that 
are also infalling and normalizable-where normalizable 
in this case means that the term x^drC'-^'* vanishes at 
infinity (as explained in the discussion around 

So far it appears that we have found a solution that is 
both infalling and normalizable, for all small frequency 
and momenta. This should not be possible, and indeed 
we have yet to impose the constraint ([70]) . We find then 
the relation 

2 2 1 2 2 '-^x /■oo\ 

uj =v,k V, = — . (83) 
at 

This is consistent only because of the radial independence 
of ax and at, which away from the boundary involves 
fluctuations of the pion field. Thus we have shown that 
there exists a normalizable infalling mode provided that 
w, k obey a linear relation. This is our main result. 



Of course as the force ttq vanishes the only solution that is in- 
falling and normalizable becomes the trivial one A = 0. 
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It is convenient to evaluate (|8T|) and ((82)) at r = oo, 
where the terms depending on C*'^ do not contribute, 
which leads to 



Qfj; = To- lim v-gg g Ora^, 
gX 

at -To- lim y^g'"^g''*drat 
gA 



(84) 
(85) 



Note that this is equivalent to the statement that at the 
boundary the right-hand side of ([70]) vanishes — in the 
absence of a Goldstonc boson source the current is con- 
served. Recall that at , ax obey the zero frequency and 
unitary gauge infalling wave equations. It was shown in 
jssi that on such a field configuration at infinity the ratio 
of the current Jx^t ~ drax^t to the field ax.t (which are 
unity in this case from ([79)) ) itself is simply the field the- 
ory Green's function for the current. We thus conclude 
that 



a, ^Ai?2c(w = 0,fc = 0), 
at = -XR^Gttiuj = 0,fc = 0) 



(86) 
(87) 



where 



Gf,^{uj, k) = (jp(w, k)j„{-uj, k)) 



retarded 



are momentum space retarded Green function of the spin 
current along a symmetry broken direction. We thus 
find that the spin velocity Vg can be written as 



Gxxiuj = o,fc = 0) 

-G'„(w = 0,fc = 0) 



(89) 



This is the expected expression for an antiferromagnet. 
In particular, recognizing that j57l - l60| 

Ps = Gxx{oJ = 0, fc = 0), x± = ~Gtt{io = 0,fc = 0) 

(90) 

where ps is the spin stiffness and x± transverse magnetic 
susceptibility, we recover the standard expression 



behavior of equation (|T7)) we see that for the solution to 
be nonsingular we need 



dr 



-9 ? 

9aX 



ax 47rr 
where T is the temperature. Introducing 

1 yf^g"g''''drax 



~9\ 



ax 



from ([75)) we find Gx satisfies the flow equation 



drGx{r) 



-gg^'^'x^ 



(92) 



(93) 



(94) 



This equation should be integrated from the horizon 
(where from ([M)) we have Gxirh) = 0) to the AdS bound- 
ary, where it is equal to the Green's function Gxx{^ = 
0,A; = 0). 

Similarly to above, for at we introduce 



Mr) 



gA^t 



(95) 



This is convenient, as at the horizon = 0. The corre- 
sponding flow equation for T-Lt is 



drUtir) 



~9\ 



-ggrrgtt 



nUr)y^g''x 



(96) 



and it should be integrated from the horizon, where the 
relevant boundary condition is Utirh) = 0, to infinity. In 
terms of ([93| and ([95)) the spin wave velocity ([89p can be 
written as 



Vs = {-Gx{oo)'Htioo))" 



(97) 



Evaluation of this requires knowledge of the scalar profile 
and can only be done numerically; some representative 
plots are shown in fig. O 



(91) 



Note that equation ([76)) has a simple boundary theory 
interpretation; oc cJ^ttq, as expected for the superfluid 
part of the current density. In our probe analysis, there is 
no normal component. Also note that our analysis above 
gives a nice correspondence between the Higgs mecha- 
nism in the bulk and the dynamics of Goldstone modes 
in the boundary theory. 



C. Evaluation of spin wave velocity 

We now evaluate Gxx and Gtt to find the spin velocity. 
It is convenient to derive a flow equation for them as 
in [53} . We start with Ux] examining the near- horizon 



FIG. 9; Spin wave velocity as a function of T/Tc for various 
values of the rescaled gauge coupling qa'- Qa is varied from 1 
(lowest curve) in unit increments to 4 (highest curve). 

It is interesting to see what happens near the phase 
transition, where we have some analytic control; here we 
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have ~ (T — Tc) from earlier results; now examining 
the flow equations we see that 'Ht(oo) remains finite in 
this limit, and is in fact given (up to corrections analytic 
in (T - T,)) by 



~9\ 



1 



(98) 



where S is the spin susceptibility in the unbroken phase. 
The second equality in (|98| follows from the U{\) analysis 
in Appendix D of |53j . On the other hand for Qx we find 
from 194]) that near the transition C/^ ~ (T — Tc). Thus 
the nonlinear term in the evolution can be neglected, 
leading to 



Wc find then near Tc 
dr 



dr ' 



(99) 



(100) 



Note that above Tc, ps = Gxxiy = 0, fc = 0) becomes 
identically zero. It would be interesting to extend this 
analysis to understand what happens to the spin wave 
velocity at zero temperature near the quantum critical 
transition. 



VI. EXTERNAL "MAGNETIC FIELDS" AND 
FORCED FERROMAGNETIC MAGNONS 

We have constructed a gravity description of the spon- 
taneous breaking of an SU{2) symmetry-analogous to 
the Nccl phase in a spin systcm-and displayed the asso- 
ciated spin waves. In this section we turn briefly to a fer- 
romagnet. A ferromagnet can be understood as a system 
undergoing spontaneous symmetry breaking in which the 
broken vacuum is charged under the unbroken generator. 
The spin waves in this case possess a quadratic dispersion 

This setup is not straightforward to realize in 
holography-the spontaneous generation of a nonzero spin 
density means that essentially one now wants the unbro- 
ken non-Abelian gauge field itself to condense, without 
supplying any external chemical potentials. We leave this 
for future work. 

In this section we consider something simpler-imagine 
applying an external magnetic field H to a sample con- 
taining spins, forcing them to align along the direction of 
the field. In our setup this corresponds to picking a di- 
rection in the SU{2) space - we will pick the 3 direction ~ 
and supplying a chemical potential fj, ^ H for the gauge 
field in that direction (in this section only ^ refers to the 
chemical potential for the non-Abelian gauge field; we 
will set the U{1) chemical potential to here, as it does 
not play an essential role). It is interesting to note that 
this model has been studied before: it is precisely the 



"normal phase" of the holographic p-wave superconduc- 
tor studied by [H^l, but our interpretation is different: in 
particular in those models a U{1) subgroup of the SU{2) 
is usually taken to be electric charge, which is quite dif- 
ferent from our approach. We will show that small fluc- 
tuations around such a configuration exhibit spin waves 
whose dispersion is not gapless, but whose dependence 
on momentum is indeed quadratic. 



A. Gravity setup 

Wc consider a generic spacetimc metric (j59|) and a non- 
trivial non-Abelian gauge field profile 

Al{r) = iia{r) a(r ^ oo) = 1 a(ro) = 0, (101) 

with all other matter fields turned off. In particular the 
scalar 0° is uncondcnscd and will play no role. Wc will 
also set the gauge coupling gA to 1 as when the scalar is 
inactive gA plays no role in the classical dynamics that 
follows. The profile a{r) solves the equation 
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--dr {^/—gg'-'g''dMr)) = 



(102) 



= -^/^9"g''drAl = const (103) 

where as above J"^' is the canonical momentum conjugate 
to and gives the expectation value of the spin density 
j"^* in the boundary theory as (j"^*) = J'^*(r = oo). For 
simplicity, we will ignore the backreaction to the back- 
ground metric, requiring that both yu and J^* are small 
when expressed in units of the temperature. Note that /i 
and J^* are related by the spin susceptibility 5 



3t\ 



(104) 



with S given by 

We will now slightly perturb this solution. We do not 
expect to find a gapless mode, as we have an "external 
magnetic field" ; however at fc = we do expect to find 
a normalizable mode with uj = p.. To see the existence 
of this mode consider performing a gauge transform on 
the background poip with infinitesimal gauge parameter 
A°: 

5 All = e"^'A'=yl^^ + duM" (105) 
It is convenient to work with the linear combinations 



A^ ^A^± iA^ 



(106) 



as these are charge eigenstates under the rotations in 
the 3 direction for which we have supplied a chemical 
potential. If we assume a time dependence of the form 
g-jwt jtqj, ^-i- j^^^ ^Yi^ resulting Af can be written as 



A+ ^A{A^-uj)e- 



(107) 



where A is an overall (complex) constant. This perturba- 
tion is generically not normalizable at the AdS boundary; 
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however if we demand that Af vanish at r — > cx3 then this 
fixes the frequency and we find 

Lo^fi, A+{t, r) = A^I{a{r) - l)e-*^* = A^o(r)e-*^* 

Thus there is a normahzable mode at w = /i with the 
specified radial profile Ato{r) = fi{a{r) — 1)^''. We will 
now turn on a small k to see how this mode evolves; as 
expected for a system with a background spin density, we 
will find that the frequency has a quadratic dependence 
on k. 



B. Dispersion relation for the ferromagnetic 
magnon 

To do this, we will need the full bulk Yang-Mills equa- 
tions: 



1 
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= 



(109) 



The presence of the background gauge potential A^ can 
be conveniently taken into account by defining a gauge- 
covariant partial time derivative dt that acts in the ± 
basis as 



(110) 



~i}j.t p — i^t^ikx 



dtA+, = {dt + iAl)At. ■ 

We give all fields a spacetime dependence e 
Here f2 = w— ^ will parametrize deviation from the uj ~ fi 
solution found above. The relevant equations of motion 
are those in (|109p for N = {r,t,x). We will work again 
with the bulk canonical momenta, defined as before: 

J+* ^ -V^F+^' = -^g'^g^^drAt, (111) 
= -^F+"' = -^gg^^g'^-'drAt . (112) 



We first examine the N — r component of (|109p . which 
can be written as 



(113) 



This equation is a bulk constraint that reduces on the 
boundary to the non-Abelian conservation of current. 
The remaining dynamical equations are 
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-.drJ 



+t 
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= 



(114) 
(115) 



Only one of the components of the field strength tensor 
is affected by the gauge field background: 

F+ = d^A+-dtA+ = ikA+ -{in{a-l)-iQ)A+ . (116) 



The existence of such a mode was observed before in [sjl in the 
context of conductivity of a p-wave superconductor. 



We now want to perturb around the uj = fi solution found 
above. We thus turn on a small fc-dependence in (|108p 
and consider a perturbation of the form 

A+ = ^H^t^^^t+rkcc ^ ^^k)+---) (117) 

where Ato{r) ^ {c({r) — 1) is the previously found profile 
in (|108p . Recall that = u ~ fi parametrizes departure 
from /i and will be the small parameter in the expansions 
that follow. A nonzero Af will also excite an which 
we take to have the form 



At ^At{r;n,k)e 



-ifit —iQ,t-\-ikx 



(118) 



The other components of the gauge fields can be consis- 
tently set to zero except for A^^, which are related to 
^t,x by complex conjugation. 

To obtain the magnon dispersion relation we will ex- 
pand the above equations to lowest order in the and k 
expansion and look for solutions that are both infalling 
(or regular) at the horizon and normahzable at the AdS 
boundary. These equations will have the same struc- 
ture as in the antiferromagnetic case above; we will find 
second-order radial equations for A^ and the correction 
Ati that are forced by the known solution Ato- 

More explicitly, plugging (fTT7| and (fTT8| into (fm]) - 
(jllSp . and expanding them in powers of and k, we find 
that 

A+(r) ^ 0{k), Ato - 0(fc2) (119) 
We thus introduce the following 



Atir) = k 1 



ax{oo) 



Atiir) = k'^atiir), 



(120) 

Here by construction A^ is normahzable, and from (jlisp 
we find that ax{r) satisfies the homogenous equation 

-La, {^g^'^g^^dra,) - fi\a - ifg^'g^^a, = 0. 



(121) 

There is a forced radial equation for the correction to An 
profile ati{r), but we will not need to solve it explicitly 
to find the dispersion, analogous to the antiferromagnetic 
case where the correction to the pion profile was not ex- 
plicitly needed. 
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dr (V^g'-^g^draa) = pi{a - l)g''g^^a, . (122) 



It is of course critical to note that an infalling and nor- 
mahzable solution to this equation can always be found. 

We now impose the constraint arising from the conser- 
vation of current (|113p . Again it is simplest to evaluate 
it at the AdS boundary first; here the right-hand side 
of (|113p vanishes (as by construction we are looking at 
a normahzable solution with ^^(oo) = 0). We find the 
dispersion relation 



~ UJ - 11 = — — ■ hm 

J6t r^oo 



-99 9 



ax{oo) 



(123) 
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This is the desired dispersion relation. 

Let us briefly understand the physical origin of the 
differences from the linear and gapless dispersion found 
in the antiferromagnetic case. The dispersion is not 
linear because the background value of J'^* is finite in 
the cj — > limit, unlike in the antiferromagnetic case 
where it is proportional to uj: this means that we are 
now balancing a term of order uj against one of order 
fc^, rather than a term of order w^. The mode is not 
gapless because at infinity the gauge-covariant derivative 
dt = dt + iA^ — > —iLo + ifjL, resulting in a shift in w. 
These considerations lead us to expect that if one were 
able to create a situation with a nonzero J"^* in the ab- 
sence of a background chemical potential - i.e. a true 
spontaneous ferromagnet - one would find precisely the 
gapless quadratic dispersion of the standard ferromag- 
netic magnon. 

Note that we can interpret the expression (|123|) in 
terms of field theory quantities. Again using expressions 
from [5^, we can rewrite the ratio of drUx to Qx in terms 



of a field theory correlator to find 

(124) 

with 

G«(..,fc) = (j + (c.,fc)j-(-a.,fc))^„ (125) 

i.e. Gxx is the field theory retarded correlator for j+ eval- 
uated at the nonzero frequency ui = ^. The prefactor of 
the quadratic dispersion is consistent with the expected 
result for a ferromagnet [6]| . 

We now evaluate the constraint (|113p at arbitrary r. 
Now the right-hand side no longer vanishes, and we find 

n = oj- fi = -fk^, (126) 

where 



J3t 



M(a-l)V^5''V*a.ati(r) 



ax{oo) 



+ ati{r) 



(127) 



It can readily checked from (|12ip and (|122p that 7 is 
independent of r, and evaluated at ?■ — > 00 we find that 
7 reduces to the expression in (|123p . Again, even though 
we did not need to explicitly solve for ati(r) to find the 
dispersion, its fluctuations are essential to make sure that 
the constraint from non-Abelian current conversation is 
upheld at all points in the bulk. 



C. Evaluation of dispersion 

We now turn to the evaluation of the dispersion 7. 
This bears a large formal similarity to the evaluation of 
conductivities studied in [s^]- We consider the following 
"transport coefficient" defined at all values of r: 



with 



ifiaxir) 



(128) 



The motivation behind this name will soon be made clear. 
Q is simply related to the boundary value of this object: 



7 



ia{r — > oo] fi) 



(129) 



where we have used (|104p . From ()12ip . we find that a 
obeys a simple radial evolution equation. 



d^a — i/i 



9rr / <7 , n2 



(130) 
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(r) 



(131) 



and at the horizon the value of a is fixed by the infalling 
boundary condition to be a{rfi; ^) = S(r/j) (this is ex- 
plained in detail in fssf). 

Before considering general fi, let us consider the limit 
/i — )■ 0, in which case drcr = 0, a becomes constant in r, 
and is in fact equal to the normal (Abelian) DC conduc- 
tivity anc of a t/(l) current on this background. Thus 
we find for the "magnon" dispersion (|123p 



-ik 



(132) 



This relation is not surprising; when taking /i to we 
are removing all non-Abelian effects and returning to the 
zero density system, and (jl32p simply describes the dif- 
fusion in the longitudinal channel of a U{1) current. The 
diffusion constant is seen to be ^ as is required by the 
Einstein relation. 

For generic fj, we must evaluate this expression numeri- 
cally. In this section alone for simplicity we work with the 
normal AdS-Schwarzschild metric in coordinates with 77 
set to 0. We do not require the RN metric here; even in 
the absence of a net U{1) charge, one can imagine polariz- 
ing thermally excited "spins" with an external magnetic 
field. We do not expect inclusion of a U{1) charge to 
qualitatively change the results. We work in units where 
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both a DC and S have been set to 1. The resuhs of the 
numerical evaluation are shown in Figures [TU] and [11] as 
a function of /i/T. 




(a) Re 7 



Imy 




FIG. 10: The real and imaginary parts of 7 are shown as a 
function of n/T. Note as — >■ 0, 7 — !> — i as required by the 
Einstein relation. 



FIG. 11: Movement of'7 in th^tomplei' plane as /i/T is varied 
from to 40. 

Perhaps the most obvious feature of these diagrams 
is the 7 has a large imaginary part, corresponding to a 
strong dissipation. This is due to the fact that wc arc 
evaluating our bulk equations as a perturbation about a 
nonzero frequency uj = fi; these waves are infalling at the 
black horizon and result in a large dissipation. We ex- 
pect that if we were able to construct a true spontaneous 
ferromagnet with /i = then our bulk equations would 
be evaluated at w = and 7 would then be both real and 
gapless. 

As mentioned earlier, this model is the normal phase of 
the holographic p-wave superconductor of [U- Accord- 
ing to their discussion, at sufficiently large y (roughly 

^ 15.9 in our units) this system becomes unstable to- 
wards a condensate of the vector fields A]:'^,. It would 



be interesting to understand further such a condensate 
- in our language, a persistent spin current - from the 
magnetic point of view taken in this paper. 



VII. DISCUSSION AND CONCLUSIONS 

In this paper we studied holographic phase transi- 
tions associated with condensation of a neutral order 
parameter, both at finite temperature and associated 
quantum critical point. We also considered the embed- 
ding of the neutral order parameter into the staggered 
magnetization for an antifcrromagnetic phase. We show 
that at the macroscopic level one recovers the expected 
features of the antifcrromagnetic phases including exis- 
tence of two gapless spin waves and their dispersion re- 
lations.^^ A similar discussion in a forced ferromagnetic 
phase reveals spin waves with quadratic dispersion rela- 
tions which again agrees with field theoretical expecta- 
tions. 

Our discussions left some loose ends which should be 
studied further. The first is the phase diagram for the 
system in the alternative quantization, for which we 
only gave a qualitative picture. The phase boundary 
for the low temperature stable condensed phase should 
be mapped out more precisely. In particular it should 
be better understood what happens to the system at 
w?R^ " ^fi" -'^^ seems possible that the various scal- 
ing exponents may not obey mean field behavior there. 
As mentioned earlier while the appearance of a special 

value m?R^ = —tt; has to do with nonlinear structure 
16 _ 

of the potential, we expect this value and the qualitative 
behavior we find may not depend on the details of the 
potential as far as the next order nonlinear term is given 
by (j)'^ (i.e. dictated by a Z2 symmetry). In particular, 
we expect similar phenomena should also appear in the 
context of holographic superconductors. Another impor- 
tant question is to work out the solution for backreacted 
spacctime geometry. While we have shown that deep in 
the IR one should find an AdS2 x with a different 
cosmological constant, working out the backreacted met- 
ric for all spacetime may reveal new qualitative features. 
Also it would be interesting to extend our analysis for 
both antifcrromagnetic and ferromagnetic spin waves to 
include backreactions or at nonlinear level [56| . 

The investigation of this paper suggests various other 
interesting directions to explore. We name a few here: 



1. The nature of the quantum critical point 

It would be very interesting to better understand the 
nature of the quantum critical point and quantum phase 



As discussed earlier these features only depend on symmetry 
breaking pattern and not details of the microscopic theory. 
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transitions discussed in this paper. On the gravity side, 
along with the example pointed out in j38| , they all have 
to do with violation of the BF bound in AdS2 region, 
whose critical behavior thus falls into the universality 
class discussed in [s^l , where a critical line (in the sense of 
a (0+ l)-dimensional CFT) ends at critical point. To the 
left of the critical point an IR scale is generated, exhibit- 
ing the BKT scaling behavior. Below this IR scale, the 
system flows to different IR fixed points depending on the 
charge q of the condensate. For a neutral order parame- 
ter the system flows to a fixed point with z = c», while 
to a (2-1- l)-dimensional CFT when q is sufficiently large. 
That such distinct physical phenomena share similar IR 
behavior is striking, yet viewed from (2 + l)-dimensional 
perspective, they involve completely different phenomena 
involving complete distinct fixed points. We have not 
yet studied in detail the critical behavior from uncon- 
densed side and spatial as well as temporal correlations 
near these critical points, which should be accessible from 
our gravity treatment. 

Despite our gravity description, a field theoretical un- 
derstanding of these critical points remains completely 
obscure. Can these transitions be understood as lying 
in the Landau paradigm? What do they suggest regard- 
ing the underlying mechanism for the condensation of 
charged and neutral order parameters in theories with 
gravity duals? Can one find explicit field theoretical 
models with similar features? We hope to return to these 
questions in the future. 



2. Competition between different orders 

From a more phenomenological point of view, more 
fields can be added to model more physics: e.g. one nat- 
ural addition is the inclusion of a superconducting order 
parameter corresponding to a bulk field if) charged under 
C^(l)charge- One could also turn on a potential y((/)°,V') 
to model the interaction between an AFM and super- 
conducting order parameter. For example, the simplest 
possibility one could consider is 

F(0^v) = &^lV'l' (133) 

Note that a negative value of h will lead to mutual en- 
hancement of both types of instabilities, while a posi- 
tive value of b will make them compete. We expect a 
rich phase structure to arise from such a construction, 
which may provide a strongly coupled example for many 
condensed matter problems in which such a competition 
exists^^. 

A different observation is that in the new IR fixed point 
for the condensed phase for the neutral scalar field, since 
R2 < R2, the IR dimension ([T]) for a (different) charged 



scalar field appears to be smaller than that in the uncon- 
densed phase. This appears to imply that the condensed 
phase or associated quantum critical point may enhance 
the superconducting instability, say making a previously 
stable mode unstable or enhancing the transition tem- 
perature to the superconducting order. A precise under- 
standing of this will require the explicit construction of 
the backreacted geometry. 



3. Coupling to fermions 

Finally, our work raises a number of concrete ques- 
tions regarding the properties of fermions on such back- 
grounds. It will be important to consider coupling a 
Fermi surface to our system In an antiferromagneti- 
cally ordered state, the nature of the Fermi surface pro- 
vides yet another means to characterize the phase. In 
particular, the Fermi surface will reconstruct as a result 
of the antiferromagnetic order, in a way that is depen- 
dent on the size of the Fermi surface in the paramagnetic 
phase. For the system we have analyzed here, in which 
a continuum limit has already been taken, it will be in- 
structive to see how to take into account such effects. 
Rclatedly, the evolution of the Fermi surface across an 
antiferromagnetic quantum critical point has emerged as 
an important characterization of the nature of the quan- 
tum critical point, and it will be interesting to address 
this issue in our system as well. 

Another issue related to fermions involves their low- 
energy spectral properties. Similarly to the scalar case 
described above, one expects that the IR scaling dimen- 
sion S of the fermion operator will be different before 
and after the condensation of the scalar x- In the holo- 
graphic formulation, 5 controls the spectral behavior of 
the fermions; 5 > 1 is somewhat similar to a standard 
Fermi liquid, whereas 6 < 1 corresponds to non-Fermi 
liquid behavior j20j . Thus we find the tantalizing possi- 
bility that the condensation of the scalar and thus the ex- 
istence of antiferromagnetic order could change the value 
of 5, perhaps driving the system between Fermi liquid 
and no n- Fermi liquid regimes. 
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Appendix A: Effect of magnetic field on IR 
conformal dimension 



In this appendix we outline the derivation of ([50]) 
demonstrating the effect of a background magnetic field 
on the IR conformal dimension of a charged scalar field. 



1. Dyonic black hole 

We will need to consider the effect of the magnetic field 
on the geometry. Turning on such an external bound- 
ary magnetic field for the U{1) current dual to the bulk 
gauge Bm in (O, the corresponding bulk geometry be- 
comes that of a dyonic black hole with both electric and 
magnetic charges, 

ds^ = gMNdx^'dx"" = !—{-fdt' + dx-') + ^^ (Al) 



with 



with 



Q2 + p2 j^^ 



(A2) 
(A3) 

ro is the horizon radius determined by the largest positive 
root of the redshift factor 



f{ro) = 0, 



(J2 I p2 

M = rl+ ^ ^ . (A4) 



The geometry (jAip describes the boundary theory at a 
finite density with the charge density p, energy density 
e. entropy density s, respectively given by 



p = 2- 



Q 



The external magnetic field H and temperature T are 



47ri?^ V 3r;5 



(A6) 



/is in (|A3|) corresponds to the chemical potential of the 
boundary system. 

It is convenient to work with dimensionless quantities 
by introducing 



Q = fir^, P = hrl 



and rescaling coordinates as 



i?2 



r ^ rro, {t,x) ^ — (t,x) 
after which equation (jAip becomes 

— ^r\-fdf + dx') + -^ — 



(A7) 



(A8) 



(A9) 



7' 



and 



377 = ^2 + /i^ 



(AlO) 



(All) 



Setting h = above one recovers the metric ^ used in 
the main text. 

We will be interested in the system at zero tempera- 
ture, for which 

Q2 + p2 ^ 3^4 ^2 _^ ^2 ^ 3 ^^-^2) 

and the near horizon region becomes AdS2 x with 
curvature radius 



i?2 



_R_ 



(A13) 



Note that in the zero temperature limit, due to conformal 
invariance of the underlying vacuum theory the physi- 
cally relevant quantity is the dimensionless ratio 



Using the second relation in (|A12[) we can then express 
bulk quantities like h (and thus /i) in terms of b, 



Vl + 1262 _ I 
2b 



(A15) 



2. Scalar operator dimension in the IR 

Now consider a scalar field in AdS4 of charge q and 
mass 771, with an action 

S = - j d^xy^ [{DMCpyn'^'cj} + m^(l)*(f\ , (A16) 

where the gauge-covariant derivative satisfies 

Dm4> = {dM - iqBM)^ . (A17) 
Note that the action (|A16[) depends on q only through 



hq 



(A18) 



which arc the effective chemical potential and effective 
magnetic field for a field of charge q. 

This problem is now similar to a Landau-level analysis 
from elementary quantum mechanics. After separation 
of variables using 



''^Y{y)X{r) 



(A19) 



we find that the equations of motion can be written as 
1 



-9 



--driV^g" drX) + i-g"u^ + m'+ 5" A') X = 
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-dlY + {v" - a2) y = 



(A20) 



Let us examine a bit more closely the n = mode, 
which is the most likely to condense, 



with 

v{y) = k + hqy, u{r) 
One then finds that 



w + 1 - - 

-9u \ \ r 



Yn{y) 



(A21) 



y+-] (A22) 



with Hn the usual Hcrmite polynomials. X{r) is a radial 
profile that satisfies the scalar wave equation with zero 
magnetic field h — 0, except that the momentum k on 
each constant-r slice has been discrctizcd into Landau 
levels: 



1 

"+2 



n = 0, 1, 



(A23) 



A similar discussion can be applied to the AdS2 region, 
where one finds that each Landau level has an effective 
mass given by 



m'^ + 2\hg\ (n + 



l\ 1 



2 / i?2 



(A24) 



The rest then follows exactly from the analysis in j20l |. 
and the IR dimension for (f> is given by 



with 



mn^2 3g +4 



(A25) 



(A26) 



6 



{6\bq\ 



yi + 1262 - 1 

' 7262 



1 

" 4 
(A27) 

where we used (jAlSP to express h and ji in terms of 
the dimensionless boundary quantity b (|A14p . Also note 
that to2^2 ^ A(A - 3). This is the result ^ in the 
main text. The critical magnetic field be can be found by 
setting the quantity inside the square root to 0, and is 



D 



1 + - 2m2i?2 



-2g2 



- 12g2 



(A28) 



where 



D= {3 + 2m^R^) 



(A29) 



is a quantity that goes to when the scalar mass is pre- 
cisely at the neutral AdS2 BF bound. 

It is interesting that the expression in (|A27p containing 

b saturates at a value -^J= as 6 oo. Thus if 

2v 3 



m^R^ + y/2,\q\ < -'- 



(A30) 



no matter how large the magnetic field is, a condensate 
cannot be prevented. This is surprising and is discussed 
further in the main text. 
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